Random Sequential Adsorption of Polymer 

Michal CieslJ] 

M. Smoluchowski Institute of Physics, 
Jagiellonian University, 30-059 Krakow, Reymonta 4, Poland. 
(Dated: January 30, 2013) 

Abstract 

Random Sequential Adsorption (RSA) of polymer, modeled as a chain of identical spheres, is 
systematically studied. In order to control precisely anisotropy and number of degrees of freedom, 
two different kinds of polymers are used. In the first one, monomers are placed along a straight 
line; whereas in the second, relative orientations of particles are random. Such polymers fill a flat 
homogeneous surface randomly. The paper focuses on maximal random coverage ratio and adsorp- 
tion kinetics dependence on polymer size, shape anisotropy and numbers of degrees of freedom. 
Obtained results were discussed and compared with other numerical experiments and theoretical 
predictions. 
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I. INTRODUCTION 



Irreversible adsorption of complex particles at solid and liquid interfaces is of a major 
significance for many fields such as medicine and material sciences as well as pharmaceu- 
tical and cosmetic industries. For example, adsorption of some proteins plays crucial role 
in blood coagulation, inflammatory response, fouling of contact lenses, plaque formation, 
ultrafiltration and membrane filtration units operation. Additionally, controlled adsorption 
is fundamental for efficient chromatographic separation and purification, gel electrophoresis, 
filtration, as well as performance of bioreactors, biosensing and immunological assays. 

Since its introduction by Feder jl| , Random Sequential Adsorption (RS A) became a well 
established method of adsorption properties modeling, especially for spherical molecules. 
On the other hand, using RSA to simulate adsorption of more complicated particles, such 
as polymers or proteins, raises a question how the universal properties of RSA changes 
when it is applied to non-spherical molecules. The question has already been answered for 
basic shapes, e.g. spheroids, spherocylinders, rectangles, needles and similar {2-6]. However, 
recent studies shows that such shapes are not sufficient for modeling adsorption of common 
proteins such as for example fibrinogen |7[. Therefore, attention of investigators has lately 



been drawn to to coarse-grained modeling of complex biomolecules and polymers |8Ml2| 



This study focuses on RSA of polymers on flat and homogeneous two dimensional collector 



surface. Similar model has been investigated by Adamczyk et al. 



131 ] : the authors, however, 



have studied adsorption on squared grid. Other works in this field used different polymers 



models e. g. 14j or assumed specific conformations of polymers and used different simulation 
methods 15[. In all of them, authors have focused on proper modeling of a specific polymer 
and therefore numerical method used for adsorption modeling was treated only as a tool. 
Here the main focus is set on properties of the tool. The main purpose of this work is to find 
out how the kinetics of RSA as well as basic characteristics of obtained adsorption monolayers 
depend on particle elongation and its number of degrees of freedom when the particle shape 
is approximated using coarse-grained approach. In our study, polymer is treated as a kind 
of a toy-model of complex molecule where both shape anisotropy and number of degrees of 
freedom are easy to control by merely changing the number of monomers. Therefore the 
model seems to be the simplest, yet universal tool for determining properties of RSA as well 
as complex particles adsorption. 



II. MODEL 



Polymer is modeled as a chain of identical, spherical monomers. In this study, two kinds 
of polymer models are used: 

i: stiff - monomers are placed along a straight line. In this model, elongation is controlled 
by number of monomers, whereas number of degrees of freedom remains constant; 

ii: flexible - monomer can rotate freely around its neighbors, however, it cannot overlap with 
other monomers. It is assumed that all monomers lays in a single plane. The number 
of degrees of freedom increases with monomer count but the increase is non-linear, 
due to excluded volume effect. 

Maximal random coverages are generated using RSA algorithm, which is based on inde- 
pendent, repeated attempts of adding polymer to a covering layer. The numerical procedure 
is carried out in the following steps: 

i: a virtual polymer is randomly created in such way that the center of each monomer is 
located on a collector; 

ii: an overlapping test is performed for previously adsorbed nearest neighbors of the virtual 
polymer. The test checks if surface-to-surface distance between monomers is not less 
than zero; 

iii: if there is no overlap the virtual polymer is irreversibly adsorbed and added to an existing 
covering layer. It's position does not change during further calculations; 

iiii: if there is an overlap the virtual polymer is removed and abandoned. 

Attempts are repeated iteratively. Their number is commonly expressed in dimensionless 
time units: 

to = N^, (1) 

where N is a number of attempts, Sp denotes surface area of a single polymer projection on 
a collector, and Sc is the collector area. Here, Sp = wnr\ for polymers build of n spherical 
monomers of radius r$. Square surface of Sc = (400 -ro) 2 area was used as the collector. The 
total number of iterations in each simulation was 10 5 £o- Analyzed polymers were build of 2 
to 20 monomers. For each polymer size and type 100 independent numerical experiments 
have been performed. 
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III. RESULTS AND DISCUSSION 

Example layers obtained from numerical simulations described in SecJTl]are presented in 
FigJU The main parameter characterizing obtained layers is coverage ratio: 




FIG. 1: Examples of layers formed by flexible (upper) and stiff (lower) polymers. The 
polymers are made of 5 monomers (left) and 20 monomers (right). 



6(t)=n P (t)^- (2) 

which, after infinite iteration of RSA algorithm, approaches the maximal random coverage 
ratio Omax = 0(t — > oo). Parameter np denotes here number of adsorbed polymers. To 
effectively measure 6 max using a finite-time computer simulations, appropriate model of 
adsorption kinetics is needed. 
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A. RSA kinetics 



Adsorption kinetics can be described analytically for two cases: low coverage limit and 
jamming limit [5 , s(] . In general, probability of adsorption depends on uncovered collector 
area described by Available Surface Function (ASF). For example, when coverage is very 
low, ASF decreases linearly with a number of adsorbed particles. When coverage increases 



16]. 



the ASF decay slows down as two or more particles can block the same space e.g. 
Therefore, for low coverage limit ASF is often approximated by: 

ASF{6) = 1 - C x d + C 2 6 2 + o(6> 2 ) (3) 

where C\ corresponds to an area blocked by single molecule and C 2 accounts overlapping 



of those areas [16|. It is worth to notice that C\ and C 2 are directly connected with vi- 
ral coefficients, which can be also calculated from Meyers diagrams, describing adsorbate 
particles in thermodynamic equilibrium. In case of stiff polymer, the particle shape can be 



approximated by a spherocylinder. Then, Ci(sc) can be analytically derived as [171 ]: 



where L is convex particle perimeter and A is its area. Though spherocylinder approxi- 
mates surface blocked by adsorbed polymer particle properly, its area is slightly larger than 
the polymer model area. Therefore, C\(sc) underestimates the real value and should be 
multiplied by a ratio of those areas: 

,-1 7r + 4(n — 1) 8 / n\ „ . ._. 
Ci = C 1{sc) « - 1 + - , for large n 5 

Parameter C 2 can be obtained only numerically e.g. 18[. In case of the flexible chain, there 
are no analitical predictions for either C\ or C 2 when n > 1. 

Parameter C\ has been numerically estimated by fitting ASF (6) defined as ([3]) for 6 < 
0.29 max . For stiff polymer, values obtained comply well with the predicted (jSJ). For flexible 
one, C\ grows with the number of monomers, which can be approximated with a power 
law (see Figj2]). As expected, stiff polymer particle blocks more area than the flexible one. 
Moreover, C\ growth with polymer size is significantly faster for a linear particle. 

The second limit of RSA kinetics (collector almost maximally filled) is of major impor- 
tance for determining maximal random coverage ratio basing on finite-time simulations. 



Uio 




FIG. 2: Dependence of blocking parameter C\ on polymer size. Dots and squares are 
simulation data, whereas lines corresponds to fits: fl5]) for stiff polymer and C\ = 4.03 • n 0,22 

for the flexible one. 



Coverage growth is then governed by Feder's law 



QQC]: 



9(t) = 9 V 



A ■ t~d 



(6) 



where A is a coefficient and d is interpreted as collector dimension 19j in case of spherical 



particles adsorption or more generally as a number of degrees of freedom [21] . Feder's 
law was confirmed for RSA of spheres in several collector dimensions 22j, including non- 
integral {23), 27)]' as we U as f° r elongated particles jf], Q. The analysis of exponent in (JSJ) 
estimated using coverage ratio growth presented in FigJH] reveals at least three things worth 
noticing. First, the value for n = 1 (spheres) is close to —0.5 as predicted theoretically 
and confirmed in earlier studies []] • Second, the value for dimer case (n = 2) is significantly 
higher, which reflects more degrees of freedom for such particles. This results differs from one 



obtained earlier 



25j. It results from less accurate approximation method used there. Third, 



the stiff polymer exponent approaches d ~ 3 for n = 3, which reflects third, orientational 
degree of freedom, and stays at this level despite further increase in monomers number. In 
case of flexible polymer, the increase is also curbed, but at higher d ~ 10 value. It suggest 
that hard core interaction between monomers limits an infinite growth of the number of 
degrees of freedom with polymer length. 
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FIG. 3: Dependence of exponent in ([6]) on polymer size. Dots and squares are simulation 
data, whereas lines corresponds to exponential fits: (—0.32 — 0.54 exp[— 1.28 • n]) for stiff 
polymer and (—0.09 — 0.65 exp[— 0.49 • n]) for the flexible one. 

B. Maximal random coverage ratio 



Maximal random coverage ratios were determined by extending obtained RSA kinetics 
into infinite time according to Feder's law ((6]) and exponents estimated in Sec fill Al Results 
are presented in FigJH The analysis of RSA of dimer shows, within error margin, the same 
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FIG. 4: Dependence of maximal random coverage ratio on polymer size. Dots and squares 
are simulation data, whereas lines correspond to power fits: 0.63 ■ ri~ 019 for stiff polymer 

and 0.73 ■ n~ 0A7 for the flexible one. 
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maximal random coverage ratio for spheres as obtained for dimers 251 ] . This, almost constant 
value of the ratio is observed in the range of n < 6, but only for flexible polymer model. 
Then, rapid decay of coverage ratio is observed. Nevertheless, values obtained for flexible 
polymer are much higher (approx. 20%) than for stiff polymer. In case of stiff polymer, 
coverage ratio decay starts at n > 3. In both the cases, the decay can be fitted using power 
function. 



C. Density autocorrelation and orientational ordering 

Density autocorrelation function G(r), also known as two-point correlation function, is 
defined here as a mean probability of finding two monomers at distance r, regardless of 
whether they come from the same polymer chain or from different particles. The G(r) plot 
is presented in Figj5j In case of spheres, the density autocorrelation function for maximal 
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FIG. 5: Density autocorrelation function. Data for n = 5 for flexible and stiff polymers 
were plotted together with data for spheres (n = 1) used here as a reference level. All 
functions were normalized to a mean density of monomers and, therefore, they oscillate 

around G{r) = 1. 



■ 1 ' 

1 


| < | i | ' | i | ' | < 

i flexible (n=5) 

stiff (n=5) 

j! spheres (n=l) 

i: 
i: 


l\ 

i \ 


!j i 

i : f- 
: i ;i 

/ ,' : ; 


'1 \ 


;\ :' 

^=^.^.TB 1 :-..,» U 


1 


i . i i i . i . i 



random coverages exhibits some universal properties, such as logarithmic singularity for 



r — > 2rn 



19 



20[ and fast superexponential decay when r ^> r 



However, even for 



relatively short polymers, those properties cannot be observed, due to periodic structure of 
particles (especially in a stiff polymer). On the other hand, at distances larger than polymer 



S 



length, almost no density correlation is observed. 

Elongated particle can form orientationally ordered structures, e.g. liquid crystals. For 
RSA on infinite collector, when particles orientations are randomly selected according to 
uniform probability distribution, there is no reason for the global orientational order to 
appear. Nevertheless, forming of local ordered domains is possible (25), H] because parallely 
aligned particles require less space (see also FigJT]). Therefore, subsequent adsorbed polymers 
will align parallel to their neighbors, especially when adsorption approaches jamming limit. 
Situation changes when adsorption on finite collectors is investigated. Adsorption conditions 
at collector boundaries may promote specific ordering. To measure it quantitatively the order 
parameter has to be defined [24]. 



N 



N < 



J2( X i COS 4> + Vi sin ' 



i=l 



(7) 



where iV is a number of molecules in a layer, [xi, yi\ is a unit vector pointing from one end of 
the polymer to the other, and denotes mean direction of all particles and can be calculated 



as m 



25]. Order parameter q is normalized so as to vanish in totally disordered layers and 



to equal 1 for perfectly ordered systems. The mean values of q for obtained coverages are 
presented in Fig[6l As expected, global orientational order increases with polymer size; 
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FIG. 6: Dependence of global orientational order parameter on polymer size. Dots and 
squares are numerical data, whereas lines are linear fits: —0.00013 + 0.0040 • n for stiff 
polymer and y = 0.0044 + 0.0017 ■ n for flexible one. 



however, its value is rather small even for the longest stiff polymer because collector area is 
relatively large. 
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To study local ordering, the following function was introduced: 



q(\r\ 



-0(£+r)] 2 ) 



(8) 



where 4>{x) is a unit vector along local orientational ordering at point x and )•) > is an 
average over particle pairs at a distance r, measured as a distance between centers of the 
closest monomers. Relation q{r) shows how the local ordering propagates in a layer. As 
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FIG. 7: Local orientational ordering propagation inside covering layer build of stiff 

polymers. 

shown in FigJTJ ordering vanishes quickly for short stiff polymers. For the longest one, it 
is significantly larger than even at distance of 10ro, which reflects a tendency for parallel 
alignment. Slight minimum around r ~ 3.5 is connected with the definition of distance 
between particles, which allows two, even perpendicular polymers, to be as close as parallel 



ones. 



IV. SUMMARY 



In the study, coarse-grain model of polymer is used to test dependence of RSA proper- 
ties on number of degrees of freedom and elongation of adsorbate particles. RSA kinetics 
at low coverage limit is sensitive to molecule shape anisotropy whereas at jamming limit 
it is governed by number of degrees of freedom, which generally is consistent with Feder's 
law predictions. Maximal random coverage ratios did not change for short (n < 6) flexi- 
ble polymers, while rapid drop was observed for stiff polymer starting at n — 3. Density 
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autocorrelations for maximally covered layers reflect the inner structure of adsorbate parti- 
cles, especially for stiff polymers. Additionally, in that case, significant local orientational 
ordering was observed, which reflects domain structure of such layer. 
This work was supported by grant MNiSW/N N204 439040. 
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